An unusually brilliant transient in the galaxy Messier 85 
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Historically, variable and transient sources have both surprised astronomers 
and provided new views of the heavens. Here we report the discovery of 
an optical transient in the outskirts of the lenticular galaxy Messier 85 in 
the Virgo Cluster. With a peak absolute R magnitude of —12 this event is 
distinctly brighter than novae, but fainter than type la supernovae (expected 
from a population of old stars in lenticular galaxies). Archival images of the 
field do not show a luminous star at that position with an upper limit of 
g ~ —4.1, so it is unlikely to be a giant eruption from a luminous blue variable 
star. Over a two month period the transient emitted radiation energy of 
almost 10^^ erg and subsequently faded in the opical sky. It is similar to, but 
more luminous at peak by a factor of 6 than, an enigmatic transient in the 
galaxy M31. A possible origin of M85 OT2006-1 is a stellar merger. If so, 
searches for similar events in nearby galaxies will not only allow to study the 
physics of hyper-Eddington sources, but also probe an important phase in the 
evolution of stellar binary systems. 



On January 7, 2006, the Lick Observatory Supernova Search (LOSS) team, in their 
daily circular, reported a source with apparent unfiltered magnitude of ~ 19.3 projected 
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2.3 kpc from the center of the lenticular (SO) galaxy Messier 85 (M85; also known as 
NGC 4382)*, a member of the Virgo cluster of galaxies; see Figure 1. For reasons stated 
below we believe that the transient lies within M85 and thus we name the source as 
M85 Optical Transient 2006-1 (or M85OT2006-1 for short). Fortuitously this field was 
observed by the Hubble Space Telescope (HST) three years earher. From the archival HST 
data we derive a pre-explosion limiting magnitude of F475W=— 4.3. 

On UT 2006 January 8 we initiated an optical photometric campaign with the auto- 
mated Palomar 60-inch telescope (see Figure 2 and Table in the electronic supplementary 
material). The light curve with a plateau of 70 days is unlike that of type la supernova. 
The plateau duration is also too short for an outburst from a Luminous Blue Variable 
(LBV; 1] Carina). 

We began a program of spectroscopic observations with the Palomar Hale and the 
Keck I telescopes (Figure 3). The Palomar spectrum obtained on UT 2006 January 8 did 
not contain any strong emission feature; the spectrum could be adequately described by a 
black body spectrum with effective temperature of approximately Tcs ~ 4600 K. Likewise 
the UT 2006 February 3 Keck spectrum was also featureless but unfortunately did not 
cover Ha. The UT 2006 February 23 and 24 Keck spectra showed a similar continuum 
but a number of emission lines were readily detected (Figure 3). Since the latter spectra 
were the deepest it is likely that the lines were seen due to better sensitivity. 

We associate the strongest line at wavelength, A ~ 6587 A and that at A ~ 4874 A 
with Ha and H/3, respectively. Accepting this identification, the mean heliocentric (peak) 
velocity of the pair is 880 ±130 kms~^. We were unable to conclusively identify the 
remaining lines but do note that the spectra of many hypergiants contain a number of 
unidentified^'^ emission lines. 

The systemic velocity^ of M85 is 729 ± 2 km s"-*^ and the velocity dispersion^ in the 
vicinity of the optical transient is 200 km s""^. The peak velocity of the Balmer lines is 
thus consistent with M85 OT2006-1 being located in M85. Accepting that M85 is the 
host galaxy (for which we adopt a distance of 15 Mpc, the standard distance to the Virgo 

*see http : //nedwww . ipac . caltech . edu 
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cluster^) the absolute i?-band magnitude of M85 OT2006-1 is —12 mag. This peak flux 
is ten times brighter than the brightest nova but (at least) ten times less luminous than 
supernova of the Type la (the sort expected in a lenticular galaxy). The narrow line width 
of the Ha line, ~ 350 ± 140kms~^ (see Figure 3), argue independently against both a 
nova and a supernova (including of type II) origin. 

The Galactic foreground extinction towards M85 is negligible, ^4^ = 0.08.^ The source 
intrinsic attenuation can be derived by comparing the observed ratio of the emission lines 
fluxes of Ha [(3.2 ± 0.2) x lO'^^ ergs"^ cm'^] and H/3 [(0.9 ± 0.1) x lO'^^ergs"^ cm-^] 
and the theoretical value of 3.05^ We estimate E{B — V)— 0.14lo;i4 which corresponds 
to an i?-band extinction of 0.401o;4q mag. This is to low to explain the unusual color and 
temperature of M85OT2006-1 with a strongly absorbed nova, supernova or LBV. 

We searched archival data from HSTe, the Spitzer Space Telescope and the Chandra 
X-ray Observatory with the view of constraining the progenitor. There is no evidence for 
a bright progenitor nor do we see tracers of massive star progenitors (see Figure 1). This 
finding (and the shorter duration) rule out that M85OT2006-1 is an LBV because LBVs 
are amongst^ the brightest stars. My < —8. Along these lines we note that M85 is a 
galaxy composed of old stars with a possible trace of a spiral arm. We conclude that the 
M85OT2006-1 likely arises from a population of stars with mass of few Mq or smaller. 

We now turn to the physical parameters of M85OT2006-1. The bolometric luminos- 
ity flux (as traced by iird^ufy] here, is the spectral flux density at frequency u) of 
M85OT2006-1 peaks at Lp ~ 2 x 10^° erg s^^ Over the flrst two months the total radi- 
ated energy is about E'ph ~ 6 x 10^^ erg. The inferred blackbody radius of the object is 
substantial, R = [Lp/(47rcTBTjj)]i/2 ~ 17(Tefr/4600 K)-^ AU. 

The closest analog to M85 OT2006-1 is M31 RV, a bright event^ (serendipitously) found 
in the bulge of Messier 31 and still lacking a satisfactory explanation. The extra-ordinary 
briUiance of M85 OT2006-1 (Figure 4) makes it doubly mysterious. The Galactic transient 
V838 Mon,^ while considerably less luminous (see Figure 4), exhibit similar plateau light 
curves and redward evolution of the broad-band spectrum. 

^Case B recombination,'' low-density limit, T = 5000 K. 
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The distinctive physical parameters (relative to novae and supernovae; see Figure 4) 
and the potential connection to a fundamental stellar process (merger) ^'^ may warrant 
coining a name. We suggest the simple name luminous red nova with the adjectives 
highlighting the principal characteristics of M85 OT2006-1. Statistics (including especially 
the nature of the host galaxies) and follow up studies would help astronomers unravel the 
origin of these enigmatic transients and also study the physics of hyper- Eddington sources. 

Received 2006 May 14; Accepted . 

1. Humphreys, R. M. & Davidson, K. The luminous blue variables: Astrophysical geysers. 
Publ. Astr. Soc. Pacific 106, 1025-1051 October 1994. 

2. Humphreys, R. M., Jones, T. J., Poloniski, E., Koppelman, M., Helton, A., McQuinn, K., 
Gehrz, R. D., Woodward, C. E., Wagner, R. M., Gordon, K., Hinz, J. k Willner, S. R 
M33's Variable A: A Hypergiant Star More Than 35 YEARS in Eruption. Astron. J. 131, 
2105-2113 April 2006. 

3. Smith, R. J., Lucey, J. R., Hudson, M. J., Schlegel, D. J. & Davics, R. L. Streaming motions 
of galaxy clusters within 12000 km s~^ - I. New spectroscopic data. Mon. Not. R. astr. Soc. 
313, 469-490 April 2000. 

4. Emsellem, E., Cappcllari, M., Peletier, R. F., McDermid, R. M., Bacon, R., Bureau, M., 
Copin, Y., Davics, R. L., Krajnovic, D., Kuntschner, H., Miller, B. W. & de Zeeuw, P. T. 
The SAURON project - HI. Integral-field absorption- line kinematics of 48 elliptical and 
lenticular galaxies. Mon. Not. R. astr. Soc. 352, 721-743 August 2004. 

5. Freedman, W. L., Madore, B. F., Gibson, B. K., Fcrraresc, L., Kelson, D. D., Sakai, S., 
Mould, J. R., Kennicutt, R. C., Ford, H. C., Graham, J. A., Huchra, J. P., Hughes, S. M. G., 
Illingworth, G. D., Macri, L. M. k, Stetson, P. B. Final Results from the Hubble Space 
Telescope Key Project to Measure the Hubble Constant. Astrophys. J. 553, 47-72 May 
2001. 

6. Schlegel, D. J., Finkbeiner, D. P. k Davis, M. Maps of Dust Infrared Emission for Use 
in Estimation of Reddening and Cosmic Microwave Background Radiation Foregrounds. 
Astrophys. J. 500, 525-553 June 1998. 



Brilliant transient in Messier 85 5 

7. Osterbrock, D. E. Astrophysics of gaseous nebulae and active galactic nuclei. Research 
supported by the University of Cahfornia, John Simon Guggenheim Memorial Foundation, 
University of Minnesota, et al. Mill Valley, CA, University Science Books, 1989, 422 p. 
(1989). 

8. Rich, R. M., Mould, J., Picard, A., Progel, J. A. Sz Davies, R. Luminous M giants in the 
bulge of M31. Astrophys. J. 341, L51-L54 June 1989. 

9. Brown, N. J., Waagen, E. O., Scovil, C, Nelson, P., Oksanen, A., Solonen, J. & Price, A. 
Peculiar variable in Monoceros. lAU Circ 7785, 1-1 January 2002. 

10. Tylenda, R. Evolution of V838 Monocerotis during and after the 2002 eruption. Astr. 
Astrophys. 436, 1009-1020 June 2005. 

11. Oke, J. B., Cohen, J. G., Carr, M., Cromer, J., Dingizian, A., Harris, F. H., Labrecque, S., 
Lucinio, R., Schaal, W., Epps, H. & Miller, J. The Keck Low- Resolution Imaging Spectrom- 
eter. Puhl. Astr. Soc. Pacific 107, 375-385 April 1995. 

12. Sivakoff, G. R., Sarazin, C. L. k, Irwin, J. A. Chandra Observations of Low-Mass X-Ray 
Binaries and Diffuse Gas in the Early-Type Galaxies NGC 4365 and NGC 4382 (M85). 
Astrophys. J. 599, 218-236 December 2003. 

13. Oke, J. B. A New Cassegrain Spectrograph for the Hale 5-meter Telescope. Bull. American 
Astron. Soc. 13, 509 March 1981. 

14. Retter, A. k, Marom, A. A model of an expanding giant that swallowed planets for the 
eruption of V838 Monocerotis. Mon. Not. R. astr. Soc. 345, L25-L28 October 2003. 

15. Arp, H. C. Novae in the Andromeda nebula. Astron. J. 61, 15-34 February 1956. 

16. Capaccioli, M., della Valle, M., Rosino, L. &; D'Onofrio, M. Properties of the nova population 
in M31. Astron. J. 97, 1622-1633 June 1989. 

17. Darnley, M. J., Bode, M. F., Kerins, E., Newsam, A. M., An, J., Baillon, P., Novati, S. C, 
Carr, B. J., Creze, M., Evans, N. W., Giraud-Heraud, Y., Gould, A., Hewett, P., Jetzer, P., 
Kaplan, J., Paulin-Henriksson, S., Smartt, S. J., Stalin, C. S. &; Tsapras, Y. Classical novae 
from the POINT- AGAPE microlensing survey of M31 - I. The nova catalogue. Mon. Not. 
R. astr. Soc. 353, 571-588 September 2004. 



6 Kulkarni et al. 



Supplementetry Information is linked to the online version of the paper at 
www.nature.com/nature. 

Acknowledgments: We thank D. Frail for discussion and constructive criticism. We 
would like to express our gratitude to astronomers who maintain the NED database at 

IPAC and the data archives of the Hubble Space Telescope, the Spitzer Space Telescope 
and the Chandra X-ray Telescope. Our work has been in part by NASA, NSF, the Sylvia 
and Jim Katz Foundation and the TABASCO Foundation. 

Author Information: Reprints and permissions information is available at 
npg.nature.com/reprintsandpermissions. The authors declare no competing financial 
interests. Correspondence and requests for materials should be adressed to S.R.K 
(srk@astro.caltech.edu). 




Figure 1. Optical images of the field around M85 OT2006-f ol^tained at two epoclis. Data 
were obtained with the Low-Resolution Imager and Spectrograph (LRIS;^^) at Keck (a) on Feb 
3, 2006, and the Advanced Camera for Surveys aboard the Hubble Space Telescope (b - F475W 
filter; c - F850LP filter) on Feb I, 2003. The event is located about 30" from the center of M85 
at a=12'^25™23.82'^ and 5=18°10'56.2" (J2000). After registering the Keck image to the HST 
image (rms of the transformation was 40 mas) we were able to place the following limits for a 
pre-cursor object (progenitor star): 26.8 mag in the F475W filter (exposure 750 s) and 24.7 mag 
in the F850LP filter (exposure time 1150 s). These limits exclude an LBV^ origin (for which 
My ~ —8 mag). Furthermore, we find no evidence for young stars (supergiants, clusters and 
HII regions). An analysis of Spitzer Space Telescope Infarcd Array Camera data obtained on 
Dec 21, 2004, result in 3cr upper limits of 25, 30, 60 and 75 /^Jy at 3.6, 4.5, 5.8 and 8.0 /xm, 
respectively. X OS'S' observed M85 two hundred and twenty times over 2000-2006. We found no 
transient at the position of M85 OT2006-1 to (roughly i?-band) magnitudes ranging from 20 to 
21. No X-ray emission was detected in a Chandra X-ray Observatory observation^^ obtained in 
June, 2002, with a flux upper limit of 2.7 x 10"^ cnt s"^ in the 0.3-10 keV band. 
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Figure 2. Temporal evolution of M85 OT2006-1. (a) observed light curve uncorrected for Galac- 
tic foreground extinction and (b) ufy including foreground extinction correction of Ay = 0.08.^ 
Data for the plots are given in Photometry Table of the electronic supplement and come from the 
following sources: Palomar 60-inch (P60; gRIz), the Large Format Camera (LFC) on the Palo- 
mar Hale 200-inch (P200; zRI), the Widefield Infrared Camera (WIRC) on P200 {JHK), LRIS 
on the Keck-I 10-m telescope (gRI), Persson's Auxilary Nasmyth Infrared Camera (PANIC) 



Brilliant transient in Messier 85 9 




4000 5000 6000 7000 8000 9000 

wavelength [A] 



Figure 3. Optical spectra of M85OT2006-1. Data were obtained with the Double Beam 
Spectrograph^^ (DBSP) at the Palomar Hale 200-inch telescope (grey line, 1800 s integration, 
Jan 8.53 UT, 2006) and Keck/LRIS (blue line, 3000 s, Feb 24.59 UT, 2006). Not strong emission 
or absorption features are seen in the (native) DBSP spectrum. Specifically we place a limit 
of 6 X 10^^^ erg s^^ cm^^ for an emission line in the vicinity of Ha. In the LRIS red channel 
spectrum the brightest emission feature is at A = 6587 A (flux of (3.2ib0.2) x 10~^^ ergs~^ cm~^) 
which we identify with redshifted Ha. The velocity of the line center is 1020 it 150kms~^ (see 
inset). On the blue side, the strongest feature is at A = 4875 A corresponding to redshifted 
(700 ± lOOkms^i) H/3 (flux of (0.9 ± 0.1) x lO"^*^ ergs~^ cm"^) _ ^j^g equivalent widths are 
10 ± lA (Ha) and 5 ± 1 A (H/3). The fuh width at half maximum (FWHM) of the Ha hue, 
after accounting for the instrumental FWHM, is 350 it 140 km s^^. In addition we detect the 
following emission lines (central wavelengths, typical uncertainty of 1 A; and fluxes, unit of 
10-16 erg s-icm-2): 4115 A (0.3 ±0.1), 6428 A (0.9 ±0.1), 6527 A (1.5 ±0.4), 8079 A (0.8 ±0.1) 
and 8106 (0.7 ± 0.1). Further LRIS spectra were obtained on UT 2006 February 3 and 23 (not 
shown here). The February 3rd LRIS spectrum did not include the Ha wavelength. For this 
spectrum, using a sliding 10 A window we were able to set a 3-cr upper limit of ^ 6 x 10"^^ erg 
cm~^ s~^ in the vicinity of H/5. 
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Figure 4. Phase space of cosmic explosive and eruptive transients. The vertical axis is the 
peak brightness in the i?-band and the horizontal axis is the duration of the event (r). Events 
are represented by circles with the color at peak magnitude coded as follows: blue {g — r ^0), 
green {1 ^ g - r ^ 0) and red {g -r ^1). M85 OT2006-1 and the M31 RV^ clearly stand out 
in this figure in the following respects: (i) brighter than novae but (ii) less luminous than most 
supernovae (especially of type la indicated with a 2a brightness band) and (Hi) distinctly red 
color when compared to sub-luminous core collapse supernovae (such as SN 1987A). Finally, the 
two events, unlike LBVs and core collapse supernvoae, do not arise in star-forming regions. For 
any reasonable progenitor mass, both events exhibit hyper-Eddington peak luminosities, similar 
to the sources V838 Mon.^^'^'^ Furthermore, both sources also characterized by low expansion 
velocity (< 1000 km s^^) and a strong redward evolution of the peak frequency. For these 
objects, r, is the "plateau" time scale. For novae, r is the time scale in which the nova fades 
by two magnitudes, t2- Filled circles show the positions of 82 novae observed^^'^^ in Messier 31 
(assuming^"^ V — R = 0.56 at peak). The brightest stars in our Galaxy are highly variable but 
these objects (marked "LBV") are clearly distinguished by long variability timescales and a high 
quiescent magnitude. The dashed line {R = —7.5 mag) is the Eddington limit for a IMq G-type 
star. 



